The Generation, Evolution and Decay of Pure Quantum Turbulence: 

Biot-Savart Simulation 
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A zero temperature superfluid is arguably the simplest system in which to study complex fluid 
dynamics, such as turbulence. We describe computer simulations of such turbulence and compare 
the results directly with recent experiments in superfluid ^He-B. We are able to follow the entire 
process of the production, evolution, and decay of quantum turbulence. We find striking agreement 
between simulation and experiment and gain new insights into the mechanisms involved. 

PACS numbers: 67.30.H-, 67.30.hb, 67.30.he, 47.37.+q 



The tangle of similar quantized vortices in a pure co- 
herent condensate provides a very simple model system 
for studying turbulence, which is more amenable to sim- 
ulation and subsequent detailed analysis than its clas- 
sical counterpart. Two recent developments inform the 
present work. First, quantum turbulence was recently 
observed in superfluid ^He. This medium turns out to be 
particularly useful as there are a range of very sensitive 
well-developed techniques which can be easily adapted 
for turbulence studies. Secondly, advances in computing 
power and techniques now allow us to make full Biot- 
Savart simulations of turbulence over relatively large vol- 
umes allowing direct comparison with experiment. 

We present a simulation of grid turbulence in super- 
fluid ^He-B in the zero temperature limit, where the nor- 
mal fluid fraction is negligible, and compare with recent 
experiments. The vorticity is initially generated as a gas 
of uniform vortex loops. We are able to follow the en- 
suing collision and recombination of these loops to form 
a turbulent tangle. After stopping the loop generation 
we can also follow the turbulence as it freely decays. The 
simulation encompasses the full range of length-scale evo- 
lution from initial loops to a developed tangle with large 
scale structure which decays back to small scales via a 
Richardson-like cascade, thus providing an ideal scenario 
for understanding various aspects of pure quantum tur- 
bulence. The simulation reproduces many of the fea- 
tures observed experimentally and provides further in- 
sights into the processes involved. 

In the low temperature limit, the few remaining ther- 
mal excitations in ^He-B comprise a highly ballistic di- 
lute gas of quasiparticles. The quantized circulation 
n = hj^vfiz around vortex cores gives a very large 
cross-section for Andreev reflection of quasiparticles, 
which allows experimental studies of vortices using well- 
developed quasiparticle detection techniques which are 
only possible in ^He-B [l|, i2|]. Crucially, Andreev reflec- 
tion has no significant effect on vortex dynamics so the 
normal fluid can be disregarded, which greatly simplifies 
the description and reduces computation times. 



The simulation is directed towards recent experiments 
on a vibrating grid in ^He-B jSKZI]. At low velocities the 
grid is found to emit a cloud of independent vortex rings 
[3] . Computer simulations show that in a transverse oscil- 
lating flow, vortex lines with fixed ends can produce rings 
[8] . When the frequency matches the resonance condition 
for exciting the fundamental Kelvin- wave mode, the line 
couples strongly to the flow and it stretches and twists 
with increasing absorbed energy. Eventually it recon- 
nects to produce a free ring with a size comparable to the 
initial line length [8] . The ring propagates away and the 
remaining pinned vortex repeats the process. The ring 
production rate increases monotonically with the velocity 
amplitude of the oscillating flow [8]. Therefore, a possible 
scenario is that many vortex lines of various lengths are 
pinned to the grid mesh used in the experiments, and 
these will produce a cloud of vortex rings. The subse- 
quent evolution of the cloud of rings depends crucially 
on the ring density. At low densities, the rings propagate 
away from the grid at their self-induced velocity. Above 
a critical density the rings no longer remain independent 
but, in a cascade, collide, reconnect and form a vortex 
tangle which disperses on much longer timescales. 

To simulate the grid experiments, we follow the dy- 
namics of vortex rings injected into a simulation "cell" 
such that the left hand side of the cell represents the 
face of the grid. The experimental grid has dimensions 
5.1 X 2.8 mm and is made up of a mesh of 10 /im square 
wires separated by 50 /im. The simulation cell is nec- 
essarily smaller owing to computing restraints, and is 
a box of cross-section A = 200/imx200/im and length 
Lo=600/im, shown in Fig. 1. In the transverse directions 
the cell has periodic boundary conditions, which effec- 
tively allows for vorticity entering the simulation volume 
from other parts of the grid. The box length is compa- 
rable to the distance between grid and detector in the 
experiments. The experimental grid is likely to emit a 
range of ring sizes, limited by the mesh size of 50 /im and 
controlled by the Kelvin- wave resonance condition at the 
grid oscillation frequency of ~1300IIz, corresponding to 
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FIG. 1: Simulation of quantum turbulence formation. Each 
frame shows the vortex configuration at the labeled time. 
Rings injected from the left quickly collide and recombine 
producing a vortex tangle which evolves on longer timescales. 



ring diameters of ~ 5 /im, confirmed by measurements of 
their decay at higher temperatures by mutual friction [9| . 
In the simulations, vortex rings of diameter D=20/im 
are injected at the left-hand side of the cell at a regu- 
lar time interval r^ but at a random position and at a 
random angle within a ^20° cone around the forward 
direction. The rings travel at a self-induced velocity of 
V = 4.6 mm/s [lO|. We assume that vortex lines always 
reconnect on colhsion, as suggested by detailed calcula- 
tions of the intersection process [11]. 

Our simulation methods [l2| follow those of Schwarz 
[105. A vortex line is represented by a string of points 
along its core. The problem is greatly simplified at T=0 
since, with no normal fluid, each point moves with the 
local superfluid velocity which is calculated by full Biot- 
Savart integration. (Although the order parameter of su- 
perfluid ^He is more complex, the equations which govern 
superfluid flow of ^He-B in low magnetic fields are iden- 
tical to those of superfluid ^He, except that the vortices 
in ^He-B have larger, more complex, cores and the cir- 
culation quantum is a little smaller [13].) While there 
is no explicit decay process in the simulations, vorticity 
is lost in two ways. First, vortices can escape through 
the open ends of the simulation cell, equivalent to the 
escape of vorticity beyond the detector range in the ex- 
periment. Second, energy is lost when transferred to 



structures smaller than the numerical space resolution 
set by the average point spacing on the filament (0.5 jiia 
in our case). The justification of this cutoff procedure is 
discussed in detail in ref. [12]. 

At a low injection rate (r^ = 5 ms), simulations confirm 
that the rings travel essentially independently. Their rel- 
atively high speed results in a rapid loss of vorticity from 
the simulation cell when the injection ceases. This corre- 
sponds to the rapid decay of the vorticity signal observed 
in the experiments at low grid velocities [3] . 

At higher ring injection rates, corresponding to higher 
grid velocities, we see very different behavior as shown in 
Fig. 1 for Ti = 1.5 ms. Here, the rings immediately start 
to collide and reconnect, establishing a vortex tangle at 
the center of the box. This corresponds to the behavior 
observed at high grid velocities in ref. [3], namely quan- 
tum turbulence which decays on much slower timescales. 
The simulations are consistent with a sharp transition 
between the two regimes as observed experimentally. 

We now turn to the transient behavior after switching 
on the ring generation. In Fig. 2 we plot the build-up of 
the vortex line density (line length per unit volume) after 
commencing the ring injection. The upper plot shows 
line-densities obtained in the simulations, averaged over 
the simulation volume. The highest curve corresponds 
to the simulation of Fig. 1. At low ring injection rates, 
the rings propagate independently and fill the simulation 
volume within the ring transit time Lq/v ~ 0.14 s. The 
steady state line density for independent vortex rings is 
easily calculated. The injected vortex-ring flux along the 
forward direction is (Ar^)"^. Assuming that the rings 
travel independently and approximately in the forward 
direction, the equilibrium line density for the rings is 
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This is shown by the horizontal lines in Fig. 2. At higher 
injection rates the steady-state ring density occurs tran- 
siently, followed by a slower increase as the tangle forms. 

The lower part of Fig. 2 shows the experimental data 
for a range of grid velocities. At the lower velocities it 
takes around a second for the full signal to develop, the 
initial rise time being limited by the mechanical response 
time of the grid resonator. At the higher velocities we 
see a further longer time evolution coinciding with the 
onset of turbulence. Given the simple assumptions made, 
the correspondence between simulation and experiment is 
striking. The line density increases as the tangle forms as 
this sluggish object acts as an effective trap for incoming 
rings further increasing the line density. The balance 
between the ring injection rate and the dissipation (see 
below) determines the equilibrium line density. 

To estimate the timescale for developing the turbulent 
tangle, we note that the longest timescales are associ- 
ated with the largest vortex structures which will have 
sizes approaching that of the turbulent region jC. Assum- 
ing a Kolmogorov spectrum, the characteristic velocity 
on a length scale C is v{C) c^ C^/^e^/^(/:/27r)^/^ where 
C 0:^ 1.6 is the Kolmogorov constant and the dissipation 
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FIG. 2: Upper figure: the mean vortex line density from tfie 
simulations, for different ring injection-time intervals, steady- 
state line densities for free rings are indicated by the horizon- 
tal lines. Lower figure: the vortex line density as measured for 
various grid velocities, the tangle forms above ^ 3.5 mms~^. 



per unit mass can be written as e = C.n^L'^ jl4(, where 
the dimensionless constant C — 0.2 is inferred experimen- 
tally from the decay data [4]. We can then estimate the 
development timescale as 
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For the simulations, using C c^ Lq = 600 /im and L c^ 
6 X 10^ m~^ gives r ~ 0.4 s, and for the experiments with 
£ ~ 1 mm and L ::^ 1 x 10^ m~^ gives r ~ 2 s, in excellent 
agreement with the results in Fig. 2. 

We may also estimate the critical vortex ring line den- 
sity for the onset of turbulence. The spread of ring speeds 
along the forward direction is ~ v{l — cosO) ~ vd^ 12. To 
collide and reconnect, two rings must approach within a 
ring diameter, thus the probability of any given ring re- 
connecting per unit time is ~ iiD'^d^ l(2ATi). Multiplying 
by the cell transit time ~ L^^jv and the number of rings 
^ Li^l(vTi) gives the total probability of a reconnection 
event somewhere in the box, p ~ {t: /2A){D0LQ/TivY . 
Once two rings have reconnected, the larger fragment 
has a much lower velocity and is thus far more likely 
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FIG. 3: Simulation of the decay of vorticity. The starting 
configuration is the last frame in Fig. 1. The frame labels 
indicate the time after cessation of ring injection. 



to collide with other rings, rapidly leading to a vortex 
tangle. We may therefore estimate the condition for the 
onset of turbulence as p ~ 1, which gives us a critical ring 
injection time and a corresponding critical line density 
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This is a lower bound since it assumes that turbulence 
sets in as soon as any two rings collide, while in fact 
more collisions may be needed. Applied to the simula- 
tion, equation [3] predicts Lc ~ 6 x 10^ m~^ whereas the 
simulations in Fig. 2 show turbulence starting between 
L = 7 X 10'^ m~^ and L = 2 x 10^ m~^. Comparison 
with experiment is harder since the critical density will 
be affected by several other factors such as the complex 
three dimensional geometry and the distribution of ring 
sizes and emission angles. However the simple expres- 
sion gives roughly the correct order of magnitude for the 
critical line density, with ~ 10° and Lq ~ 1 mm. 

We now discuss the decay of the turbulence. Fig. 3 
shows a continuation of the simulation of Fig. 1 after 
stopping the ring injection. By 0.1 s the loss of the incom- 
ing rings leaves a clear volume to the left with a slowly 
decaying tangle in the center. There is a small loss of vor- 
ticity through the ends of the simulation cell but the ma- 
jority of the loss is into features smaller than the compu- 
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FIG. 4: Decay of the vortex line density. Circles show re- 
sults from the simulation in Fig. 3. Lines show experimental 
data from ref. [4]. Dashed line shows the expected late-time 
limiting behavior for a classical- like Richardson cascade. 



tational space resolution, which effectively gives a length 
scale below which vorticity is dissipated as assumed in 
the Richardson cascade of classical turbulence. 

In Fig. 4, we compare directly the decay of the line den- 
sity from the simulation of figure 3 with the experimental 
measurements at high grid velocities [4]. The late-time 
limiting behavior, L ~ t~^/^, shown by the dashed line 
in the figure, provides indirect evidence of a Richard- 
son cascade [^], transferring energy from large to small 
length scales where dissipation ultimately takes place. 
Similar behavior has been observed in superfiuid ^He ex- 
periments [15, 16] and in earlier simulations [17]. While 
in classical turbulence viscosity provides the decay mech- 
anism, turbulent decay in a superfiuid at T=0 app ears 
to be governed by the circulation quantum [% |l5|, fig , 
possibly from radiation by high-frequency Kelvin waves 
emitted by kinks left by vortex reconnections [IJ, ^M, llSl • 
At late times the simulated line density runs into prob- 
lems of the limited simulation volume, however there is 
a remarkably close similarity to the experimental data. 



Given that the simulations simply remove energy from 
small length scales with no overt physical mechanism, 
the result strongly suggests that, provided dissipation 
occurs only at small length scales, the late time decay 
is insensitive to the precise dissipation mechanism. 

In conclusion we have demonstrated that for a pure 
quantum fluid carrying turbulence we now have the com- 
putational tools to undertake a full Biot-Savart simula- 
tion of the evolution of the vorticity over significant vol- 
umes and time scales to make meaningful comparisons 
with experiments. The agreement with the observed be- 
havior in the present case is quite startling. The partic- 
ular experiment chosen for this comparison, grid turbu- 
lence, exhibits not only a Richardson- type cascade from 
larger to smaller scale structures but also a reverse cas- 
cade from smaller to larger length scales as the initial gas 
of small precursor vortex loops recombine to form larger 
and larger structures. The simulations allow a detailed 
investigation of the various processes involved. Further- 
more, the experimental observation of a vortex-ring gas 
as a precursor to turbulence was completely unexpected 
but is easily understood in the present simulations. In 
this paper we are limited to static figures. Many more 
details, such as the influence of Kelvin waves produced by 
recombination events, are revealed by studying the full 
ti me evolution (the full simulations can be seen as v ideos 
at http://matter.sci.osaka-cu.ac.jp/^bsr/tsubotag'). We 
expect that many further interesting results will arise 
from more detailed analysis. For instance, we have re- 
cently been studying fluctuations in the vortex line den- 
sity which show a Kolmogorov-like, —5/3 power law, fre- 
quency spectrum as found experimentally [6|, |20| . Details 
will be reported elsewhere. 

The ability to apply realistic simulations to experi- 
ments represents a significant milestone in the study of 
quantum turbulence which will greatly enhance our inter- 
pretation of experiments and our detailed understanding. 
We are also learning that quantum turbulence carries a 
strong resemblance to classical turbulence and may there- 
fore provide new insights into turbulence in general. 
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